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ABSTRACT 

An average of all high resolution UV spectra of supernova 1987A recorded by lUE a few days after the 
explosion reveals numerous interstellar absorption lines between 1200 and 3100 A. Features produced in 
the neutral medium of our galaxy and the LMC (e.g. those arising from O I, C I, C II, Mg I, Mg II, Fe II, 
Mali, Si n, Ni II, etc.), together with lines from much higher stages of ionization (A1 in. Si IV and C 
IV), were recorded with far greater fidelity than for previous observations toward stars in the LMC, owing 
to the supernova's extraordinary brightness and very smooth continuum. 

Although profiles are shown for a large number of species, our study emphasizes the absorption 
produced by the more highly ionized species including A1 III, Si TV and C TV. The observed absorption 
line velocity structure which covers the heliocentric velocity range from -30 to 330 km s' 1 is very 
complex. In our discussions, we assume that the absorption with v < 120 km s ^ occurs in the Milky 
Way and its halo while the absorption with v > 190 km s’l occurs in the LMC. The line profiles for A1 
ni. Si IV and C IV are converted into plots of optical depth (and column density) versus velocity. In the 
velocity range which we believe pertains to the Milky Way , we find that the profiles for Si TV and C IV 
are quite similar and much smoother than the jagged appearing A1 HI profile. On relating column 
densities, we find that while the C TV to Si TV ratio is relatively constant (ranging from 3 to 5) over the 
velocity range from 0 to 100 km s’l, the C IV to A1 m and Si IV to A1 ID ratios vary by nearly a factor of 
10. This suggests that the C TV and Si TV along this sight line in the Milky Way and its halo may have a 
common origin which differs from that for A1 HI. A1 III is almost certainly produced by ultraviolet 
photoionization in the gas of the Milky Way disk and low halo. A different origin is therefore suggested 
for Si IV and CIV. 

Absorption lines created by gas in the immediate suixoundings of the supernova are hard to differentiate 
from the very strong features created by foreground, more generally distributed material in the LMC. The 
best line for highlighting gas very near the supernova, the 1264 A line fix>m Si II in an excited 
fine-structure level, indicated that out to about 10 pc from the supernova the density of a mass-loss wind 
from the progenitor in its red giant phase is consistent with a total flux divided by the wind velocity M /v^ 
<1x10*^ M 0 yrl (km s'^)"^ . With the small accumulation of gas from such a wind, it does not appear 
possible that the pulse of energetic photons created during the supemovas's shock breakout could explain 
the strong Si TV and C IV features at LMC velocities. 
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I. introduction 

^reoccurrence of the brigh.es, supenrovain 384 ,ears (SN .987A) has provided a unique 

„ppor.„ni.yfor,hes.udvofin.ers.e,.argasbe.ween.heSunand.hel^geMagei,anic^ 

,,MC,.The50hpcsigh.Une.o.heLMCsan.p.esin.ers,enarn,a,.erwhic^ 

.a,.shgas,.o.dens.h...Wa,ha,ogas.o.densi.^^^^^^ 

andgasin,hecircun,s,eUarenvironmen,of,hesupemovauself.Thefavorab,e 

,MC ofabou.280iuns-> Dopp.er shifts absorpUon features associa.ed^U. .he LM ^ 
.Hosebe,onging.o.heM.h.yWay. Per a relatively shor. period after i.s drsorver,, SN 98 

was veo- bright in the far- and middle-ulftaviole. wavelengths accesstble to the specnogr 
aboard *e In.emanonal Uteaviole. Explorer dUE) Satellite , Through the efforts of e snt a 

USandEurop.au lUEcon™ls.tions,theh.ghdispe.ionsp.cfta,istedinTablelwereo^ 

Cl.nUa.sisof.heseda.aCdeBoere..an9S, 

;lL.„gly:hcharacterof.e — .edtutn^twee^ 

established through earlier lUE observations of bright EMC starsby Savage and de Boer ( 

1981). AnatlasoflUE in.erstehar line profdes for the sigh, line » SN1987A basedon e 

European data is found in Blades et al. (1988). 

■rre optical interstellar absonrtion line da. now available for SN 1987A are antong 
available for any inters.ehar line of sight and .veal the weald, of infomtanon whtch potent, y 
resides in such da. (Vidal-Majar et al. 1987; Magain 1987; Penini e. al. 1988) Whde e 

apecnahaveonly.odera.epho.otnedicquah.andsubs.andafty lower resolun^^ 

spectra nd y mveraec is that it includes a very 

observations from the ground, a principal advantage o 

1 a nnmhpr of the more abunda.nt elements. 

wide range of ioniration and encompasses a arg 

supernova is a pardcularly good baclcground soume, since it pmvides a relattvely smooth 

Zpgnaundcondnuumagains.wh.chnarrowa.somedmesbroad — 

“L, 

As a result, un measures of absorption 

much smaller for large sections in the supernova spectra than typ 

C de Boer (1981) or in the more recent work of Savage and M 

toward LMC stars by Savage and de Boer (ly ) 

y In Table 1 in order to obtain UV absorption 
In this paper we have averaged the spectra hs 

^ ^ ^ ■ ,^fc\nQ87A All of these spectra were 

line profiles of the highest possible quality in the duecuon 



recorded before d.e supernova had faded to dre point where conhibudons front the adjacent early 
type sta« becan* important ( Blades et al. 1988). nte details of me averaging process and plots of 
me msulunt profiles are found in 5... Tme analysis discussed in |in and 5W is mosdy concerned 
wim me new infonnation pntvided by mese data about the origin of me intersteUar absorpnon hues 
produced by C IV. Si IV and Al in at bom Milky Way and LMC velocides. 

■ntrough the analysis of an absorpnon line which is unlikely to have snong contributtons from 
tomground material in me LMC, our observations can place an upper limit on me amount of 
circumstellar ejecta from the supernova's progenitor. This Umit, in mm, indicates mat very hide of 
the strong Si IV and C IV absorption is created by a flash photoioniaation of nearby matenal at e 

time of the supernova's shock breakout. 



II. OBSERVATIONS AND REDUCTIONS 

The interstellar absorption line observations presented here were obtained with the two echelle 
spectrographs aboard the lUE satellite. The spectrographs , their respective SEC detectors and 
UV to visual light converters cover the far-UV (1 150 to 2000A , SWP camera) and middle-UV 
C1800-3300A, LWP camera) regions of the spectrum with spectral resolutions of approximately 20 
to 25 km S' 1 (FWHM). For details of the lUE satellite and its spectrographs see Boggess et al. 
(1978a,b). The individual spectra were processed at the Goddard Space Flight Center or at the ESA 
lUE analysis facility in the standard way with version 2 of the lUE software (Tumrose and 
Thompson 1984). 

Table 1 lists the spectra for which multiple high dispersion lUE spectra have been averaged in 
order to obtain interstellar line profiles with higher signal to noise than is available in single spectra. 
The various spectra were all obtained with the supernova centered in the large (10"x20") entrance 
aperture. As a result, the alignment of true spectroscopic features with respect to detector fixed 
pattern noise structure will be approximately the same for all the exposures with a given detector. 
Therefore, the averaging process will not appreciably reduce the effects of detector fixed pattern 
noise. Table 1 lists the exposure time and the date and time of the observations which were obtained 
at both the NASA and ESA ground stations. The supernova faded very rapidly during its early 
evolution particularly at far-UV wavelengths. Therefore, the various spectra differ in quality and in 
wavelength coverage. The comments to Table 1 indicate spectral regions with weak signal levels 
and also spectral regions with detector saturation problems. In producing average spectra we have 
rejected data in saturated spectral regions and in spectral regions with very low signal levels. Most 
of the final profiles were determined by averaging the data from more than two images as indicated 
in the column of Table 1 that lists the wavelength regions of the various spectra which were used in 
the averaging. 

The various spectra were averaged by first establishing the velocity shift required to bring the 
spectra into a common velocity system. This was done with reference to the many interstellar lines 
seen in the spectra. The zero point of the velocity scale was taken to be that provided by images 
SWP 30379,30381 and 30383 at X < 1800 A . The velocity shifts of the long wavelength data were 
obtained by making the middle UV lines of Fe II have the same velocity as the Fe II XI 608.456 A 
line. The velocity shifts actually applied to the individual spectra are listed in Table 1. Velocity 



shiftsofupto-18kms-l were required to register the LWP images to the average of SWP 

30379.30381 and 30383. Within a single image there was no attempt to allow for different velocity 
shifts from line to line with the exception of the lines of Si II J.1808 and A1 HI U1854 and 1862. 
Because of detector saturation, the data for these lines were obtained from a number of images 
which had weak signals for X <1800 A. The registration of the data for these lines was to the scale 
estabUshed by SWP 30379 . We note through our efforts to bring these data into velocity 
registration that 2o velocity uncertainties of about 8 km s'l exist in the combined data . Differences 
this large occasionally exist between our presentation of the results and the atlas of Blades et al. 

(1988) which are based on the European data alone and for which no special corrections were 

applied. 

In the averaging process the measured signal (i.e. flux x integration time) for the individual 
spectra are linearly interpolated into 5 km s'' velocity bins and added. The wavelength 
uncertainties introduced by the averaging process degrades the spectral resolution of the final data. 
We estimate that the resolution of the averaged data is approximately 25 km s'l (FWHM). 

For selected lines of particular interest (e.g. Si II* X 1264.730 A) where even the averaged data 
are vety noisy, we have smoothed the averaged spectra by applying a triangular weighting function 
to three adjacent 5 km s'* velocity bins with a weighting funcnon of 0.5.1. 0 and 0.5. 

In the processing of the individual specua. the linearized data numbers extracted from the 
image were corrected for the instrumental background, the echelle blaze functton and the 
instnimental sensitivity using the standard proceedures described in Tumrose and Thompson 
(1984). For absorption lines with X < 1 300 A for the SWP spectra and for X < 2300 A with the 

LWR spectra, where the crowding of echeUe orders causes this ptoceedure to over estimate the 

appropriate background . an additional level shift was determined with reference to the cores of the 
very strong interstellar lines in these wavelength regions. Ibis additional shift typically amounted to 
10 to 20 percent of the continuum level. 

The velocity structure of the intersteUar absorption to the LMC is exceedingly complex and not 
easily summarized with words or tables. We have themfore provided plots of the average observed 
flux from the supernova in the vicinity of many of the interesting intersteUar lines which are within 
the detection capabiUties of the lUE. Tbe lines plotted and echeUe orders from which the data are 
obtained are listed in Table 2. In some cases, more than one absoAing species produces a line in 



each short region of the spectrum plotted as indicated in Table 2 . The line profile plots on a 
heliocentric velocity scale are shown in Figure 1 . For the direction of the supernova, one should 
subtract 16 km s'l from a heliocentric velocity to obtain the velocity in the local standard of rest. 
Detector reseaux are identified on the profile plots with the letter R. 
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m. RESULTS 

The UV interstellar absorption line spectrum seen toward SN 1987A contains information about a 
large number of phases of the interstellar medium in both the Milky Way and in the LMC. These 
phases include those producing such neutral gas features as O I , C I, Mg I, and Cl I to the highly 
ionized gas absorption phases responsible for the lines of A1 HI, Si IV , C IV and N V. In addition, 
the absorption features of C n. Si n and Mg D with their large f values are sensitive to intervening 
gas having extremely low column densities. The emphasis in this paper is on the properties and 
origin of the absorption features produced by the most highly ionized gas found at both Milky Way 
and LMC velocities. 

a)Velocities 

The complex structure exhibited in Figure 1 appears to be a general characteristic of the LMC 
sight line, which passes through a large number of different absorbing regions in the Milky Way and 
in the LMC ( de Boer and Savage 1980 ; Savage and de Boer 1981; Songaila et al. 1986). 
Vidal-Madjar et al.(1987) noted in their very high signal to noise and high resolution (FWHM = 3 km 
s-1 ) spectra the presence of possibly 23 distinct absorption components in the heliocentric velocity 
range of 0 to 300 km s'^. Some of these many components may represent sub-structure within main 

components. 

The principal components seen in the lUE data for low ionization species occur at heliocentric 
velocities of approximately 15, 65, 120, 170, 220 and 280 km s'l . These velocities actually refer to 
features which are blends of components, a result which is apparent with reference to the higher 
resolution Ca II line data. For example the 1 5 km s' 1 component seen in the lUE data for low 
ionization species such as Si B , Mg I, and Fe H is actually probably composed of components near 
9 ,16 and 23 km s'^ while the 65 km s'^ feature is probably composed of absorption components 
near 40, 57, 65, 72, and 75 km s'l and the very strong feature near 280 km s‘^ probably contains 
seven subcomponents spread over the velocity range 250 to 295 km s'l with the strongest 
components near 280 and 286 km s'l (Vidal-Madjar et al. 1987). The velocities infeixed from the 
lUE data exhibit differences from one ion to the next. Some of these differences undoubtedly are 
produced by ion to ion differences in the relative strengths of the many components created by 
differences in the excitation conditions, elemental abundances (and their gas phase depletions), and 
perhaps the degree of line saturation. 



The low ionization lines which have large f values, such as those of C II, Mg II and Si II, are 
extremely strong and wide. For example, the Mg II X 2795.528 line exhibits nearly total absorption 
between -10 and+300 km s‘l and has half intensity absorption velocities, v. and v+ , which range 
from -30 to -1-330 km s' V It is apparent from these strong lines that there exist low ionization 
components of absorption near approximately -15 and -t-315 km s'^ which are not seen in Ca II 
absorption. 

The higher ionization lines of A1 m. Si IV and C IV exhibit absorption with components having 
their peak absoiption centered near heliocentric velocities of 20, 65, 200 and 280 km s'^ . The very 
strong feature centered near 280 km s'^ is certainly multicomponent in character with two or more 
components absorbing in the range from 250 to 300 km s'l being required to explain the breadth and 
shape of the feature. While weaker absorption occurs at other velocities, its component structure is 
not well enough defined to make velocity assignments possible. In the velocity range from -50 to 
-1-200 km s'l, the high ionization C IV lines have a distinctly "smoother" appearance than do the lines 
of A1 in and Si IV or the lines of species in lower ionization states having roughly similar absorption 
line depths (e.g. Si II ?il808.00. Mg I X2852.127 , Fe H X1608.456). A more quantitive 
comparison of the similarities and differences between the absorption profiles for the highly ionized 
species and some of the lower ionization species is found in §nib. 

The lUE data in the region of the N V doublet at XX 1238.821 and 1242.804 are extremely 
noisy. Useful information on the presence or absence of interstellar N V can only be obtained by 
processing the individual lUE echelle orders containing the lines of the N V doublet with a spectral 
extraction routine that properly corrects for the detector fixed pattern noise. Until this work is done, 
we consider any statements in the literature about the presence or absence of interstellar N V 
absorption toward the supernova based on the spectra we have examined as quite tentative. 

b) Profile Analysis Technique 

In a conventional analysis of interstellar absoiption lines, one identifies specific velocity 
components, measures their equivalent widths, and then derives column densities using a curve of 
growth based on the saturation behavior of a simple Gaussian distribution of radial velocities. 
However, for the features seen in the supernova spectrum, the principal, seemingly distinct 
components are sometimes blended together by the instrumental profile. Also, as we pointed out 
earlier we know from high resolution recordings of the absorption features from K I, Na I and Ca II 
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in „e s„pen,ova spec^m a. visible wavelengths hta. single peaks in the lUE speeha are achtally 

composed of complex, narrower components. 

one nhgh, be tempted to tty to use dte high resoludon data to create a detailed model that 

reconstructs how the ultraviolet lines should saturate as they pow sdonger. Except tn specal 
chcumstances. this is a hazardous undertaking for several reasons. First, dtfferent atomtc spec.es 
exhibit different patterns of absorption. The Ca II lines toward the supernova are more numerous an 
somewhat less "spiky” than the Na I or K 1 lines. This is probably a consequence of drastic 
reductions in the calcium depletion for components which have been shocked and accelerat - 

phenomenonrecognizedlongagoformatetial Within titegalacticplaneCRoutiyt^ 

second, the fractions of atoms existing in stages of ionization below titose favored for reg.ons 
scale with tire local election densities (and inver^ly witi. tite density of ionizing radiation). Because 

titeinterstellar medium is heterogeneous, the curvesofgmwtirofsuchUnescandti^^^ 

atoms in titeir favored (and higher) ionization stages (Spitzer and Jenkms 1975). F.nally, 

cur^ryexarrunation of the lUEdatafor ions with differingexcitationpotentials prov.de d^^ 

evidence that no single velocity pattern is globally appropnate. 

Forthereasons given above. wechosetopresentandinterprettitelUBda. in the fouo™^^^ 

graphical fashion. After defining a continuum level 1,(X) . we convened restdual tmensmes 

into apparent optical depths according to the relation 

x(X) = ln[Ica)A(X)]. 

„ ts imporumt to note, however, titat the r(k) defined above conectly depicts a smoothed optical 

deptitonly under tite conditions titat x(k) « 1 or else the velocity structure is fully resolved by e 

instrument. . 

AS a start, we shall examine the behavior of the velocity structure exhibited by an ton « 

exp^tedtobeabundantinordinaryHlandHlIregions. Rgure 2 shows loganthmtc plots o e 

apparent t 's of the column density based on the three absorption lines of S. II. In the plot, t e curv 
for different lines are displaced verticaBy on the logarithmic scale so titat each would deptc. the 

logarithm of a column density per unit velocity 

log [ N ( V ) 1 = log I - log f + 14.576 (cm -2 (km s-l) 1, 

,.betis were indeed valid. Tbefvalues for thetransitions are those givenbyShuti 

,pg,.Thefacttha.thecurvesforthestionglinesatl3()4andl526Aatesignific^^^ 
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>, 1- «1808Atadicatesthatunresolved.saturatedprontestiucturesarepresent. 

the weaker me , i j hehave if lUE were observing absorption 

. » cJrmilated how the X plots would behave it ivc. w 

intensities to various powers, snt ^ abundance front 

,eco„ve.ionsaccordtn::^— 

one vcloctty complex ^ log fX, This eoncordanee indicates that the 

those shown in Fig. 2 for ide H„nlicated by Si n. If we 

11 complex and ragged character of the Ca H velocity profiles is duph 
generally complex an gg , ^ v „f Si n we also learn from our 

accepttheCanprofilestructuresasavalidmode or os _ 

sitnnlationthatforv<220kms-l the log (r ) depicted by U-e weak 

representation for the true logC N / a ) (i.e., « ts not just a lower Ittrat). 

Tbe Shoulder locatedat.40 km s-liurhelSOSAlineseemstohavenocountetpan^ 

• • of s T at 1807 31 1 A , its apparent velocity with 

j«6AUnes. , 3 „„„,absorpdon by neutral atoms at about .75 

respect to the Si n line should be -116 km ■ i locyi 

respect xo me . e for T or Na I (Vidal-Madjar et al. 1987). 

km s-1 is not indicated by the high resolution data for K 

c) Interpretation ofM HI, Si IV and CIV Profiles 

WenowfocusontheproftleshapesinthellBrecordingsofabsontUonfea.^ 

highlylo^aedspecles. Eigures3to5showtheophcaldepUtpl«sforAl.n.S^ 

respectively. For each ion, the nansidon probability of one line is nvice that o, the 

T e ofthe weak Alin 1862.795 A line up by 0 . 3 . in order to compensate for deferences 

moving log t of the weaK velocity 

reconstruction of the pron , our data reduction. 

.n hv the addidonal smoothing and bmrong in our data reouc 

calculated the modificanon caused y 


•»= high resolu,ionda,a for Can, differences are s™a„er .han , ho 

Ihe actual optical depths are considerably less. . '' ^ 

Between -20 and 100 km c-1 tKr. 

"o.hcarlyas„elldennedinSu;::~^ 

*e inferred log N(v) between the tr in 

S iN^v; between the strong and weak Unes for Si IV or P TV . u , 

strong. These results indicate that the i • • ®'^en though the lines are 

co«— 

overveloci.shouldyieldcoluntndensideswl'rr:;!;" t 

tfte sntoothing of untesolved, saturated smtctu^s Table 3 lis T^h '' 

o--.n.is™a„nerorert.eerel.i,.„..::;:^^^^^^ 

and(3) v> 190kms-l Th,* i ’ (2) 120< v < 190 km s'^ 

• The lower and upper limits ("1.1.” and "n 1 " 

reptesentintegtattonsalongprolllesconnecdngtheloweroruppersetsofsntall ' 1 ““““ 

Enors in the condnuun, level were assented to be equal to the “ 

aboutthecondnuum,dividedby(p/3)l/2 r . ™dev.at,ons of intensity 

Ihe height and shape of the ' "tensity points used to define 

enot.) w ‘^'^"■“^^"“““■‘-“fP^-hecoherencelengthofthe 

errors.) Weestunated that the loenor in the backero„„d nn«"ofthe 

me background corresponded to about of i , 
continuum level Thf* ],tv,Wo ■ l . lu aoout o /c of the local 

teladvely strong (intervals 1 and 3) these etrors ar 1 

--groundandcondnuu.plac;::”^^^^ 

~..endd.e.atwecanon.yes.ate:r:::::r""“^^ 

and mtLXlTcrr d“^^ 

N(Si 1V,/N(A1 ,n>, N(C IV, W ZZTZi “'f 

Rgmesd 7and8 Wh-t 1 , of velocity are shown in 

son a 0 . 7 and 8. While the contparisions of Si IV and C IV with At m , h , 

is interesting how dte tado of Si IV to C tv • , • " “ 

no ot St IV to C IV ts relatively constant with velocity 
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IV discussion 

various absorption components on ,Vc 

discusses LMC gas, and §Wd discusses t accessible to the 

these discussions wUi be on the lesuits le anng ^„„stellar environment we 

,UE (e.g. A1 m. Si TV and C IV). However, rn the case 

also consider the excited fine structure levels of St n. 

a) Physical Location of the Absorbing 

Becauseoftheldnetnadcalcotuple..^ 

is difficult to unatnbigously reasonable to assign components 

absorption components detect me s'l to absorption by relatively 

near v„ah„ of about 15 km ^ 280 tan s'l wtth 

local MilkyWaygasandthecomponentsneartheLMC ^.o^ediate 
. V. TToir The uncertainty concerns the componeni 

abstuptionbygasindreL . j 2 o , nO . 200,220 km s->. In their studies of s.x 

heliocentric velocities of approximate y . j9gl) and de Boer and Savage 

stars in the direction of the LMC, Savage and e , _ „oar 65 and 120 km s'l and 
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N V toward the SN are too noisy to be of any value and lUE is unable to operate at those 
wavelengths needed to record O VI. Therefore, our assessment of possible changes in the behavior 
of line profiles will be restricted to A1 HI, Si IV and C IV. 

In §nic we noted that between -20 and 120 km s'^ the prominent double peaked absorption 
structure seen for Si II and A1 HI is not as well defined in the Si IV and C IV profiles. This change 
appears to be the result of a progressive increase in the overall breadth of the features from A1 m to 
Si IV to C IV. These differences are too large to be caused by changes in the Doppler broadening 
among the three ions at their favored temperatures from collisional ionization. For instance, the 
peak equilibrium abundances of C IV and Si IV occur near T = 10^ and 0.8x10^ K, respectively 
(Shapiro and Moore 1976) . At these temperatures, thermal Doppler broadening of C IV and Si IV 
would produce gaussian broadened line profiles with FWHM of 20 and 1 1 km s' respectively. 
However, in non-equilibrium cooling as expected in a galactic fountain flow, the recombination of 
the gas proceeds more slowly than its cooling and the temperatures at which Si IV and C IV peak in 
abundance are somewhat reduced. Clearly a broadening of less than 20 km s'^ for C IV and less 
than 1 1 km S' 1 for Si IV would be inadequate for explaining al the observed profile differences 
between Al HI, Si IV and C IV. 

The most surprising effect seen in the detailed profiles shown in Figures 3 , 4 and 5 is that the 
Si IV profiles more closely resemble the smoother appearing C FV profiles than the sharper Al HI 
profiles. From the theoretical considerations referred to above, one might expect the profile of the 
two species whose origin is expected to be dominated by photoionization (e.g. Al III and Si IV) to 
be roughly similar and to differ substantially from the profile for an ion expected to have an 
important contribution from cooling fountain gas (e.g. C IV). Additional concern for the behavior 
of Si IV with respect to C IV comes from an investigation of Figure 8 which shows the ratio of 
column densities as a function of velocity. In the velocity range of strong Milky Way absorption 
(i.e. from 0 to 100 km s'l) the ratios of C IV and Si IV column densities range from about 3 to 
about 5. If the physical processes producing the C IV (e.g. cooling fountain gas plus some 
photoionization) differ substantially from the processes which produce the Si IV (e.g. almost 
entirely photoionization) it is difficult to understand why the ratios of these ions should be so 
constant with velocity. While it appears that all the important atomic physics is included in the 
calculation of expected column densities in cooling fountain gas by Edgar and Chevalier (1986), we 
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order to begin to answer these questions we must first consider what is known about the gaseous 
environment of the SN and what absorption features might be expected for such an environment 
before the SN explosion occurred. 

A visual inspection of deep, nairow-band images of the 30 Doradus region of the LMC 
obtained in the lines of Ha and [O III] X5007 by Lasker (1971), Elliot et al.(1977) and Lortet and 
Testor (1984) reveals that SN 1987A occurred within the boundaries of the southwestern extension 
of the 30 Doradus nebula. The gaseous emission at the position of the SN is weaker than at the 
center of the nebula which lies 21 arc minutes away. However , the emission is still quite strong 
and exhibits a complex appearance which is a characteristic of the gas in 30 Doradus. An excellent 
CCD image of the [O m] emission immediately around the SN obtained using a coronograph is 
shown by d’Odorico (1987) who comments that the precursor star appears to have been located 
near the edges on a bright ridge of [O HI] emission which lies at the boundaries of a bubble. 
Although nebular emission line spectra have not yet been published, Danziger et. al. (1987) noted 
that the emission near the SN has the characteristics of a low density H II region with weak nitrogen 
lines. 

The presence of [0 10] emission reveals the highly ionized character of the gas in the general 
vicinity of the SN. The creation of O III requires 35.1 eV which is similar to the 33.5 eV needed to 
photoionize Si HI to Si IV. This implies that one might expect to measure substantial column 
densities of Si IV and perhaps even C IV along the sight line to the SN as a result of absorption by 
pre-existing H II region gas, depending, of course, on the exact position of the SN with respect to 
the dense clumps of H II region gas and whether or not the SN is in front of or behind these 
clumps. Information about the line of sight position of the SN can be inferred from the reddening 
to the SN which has been estimated to be E(B-V) = 0.15 (Fitzpatrick 1988) . Approximately 0.07 
magnitudes of this reddening is produced by the Milky Way ( McNamara and Feltz 1980) with the 
remaining 0.08 magnitudes produced in the LMC. With the LMC hydrogen to E(B-V) ratio 
appropriate for the 30 Doradus region from Fitzpatrick (1986) of 2.4 x 1Q22 atoms cm'^mag'^ , an 
LMC reddening E(B-V) = 0.08 implies the SN is beyond about 1.9 xlO^l atoms cm*2 of H I. The 
H 1 21 cm emission profiles of Rohlfs et. al. (1985) for positions near the SN give a total LMC H I 
column density of 2x1 0^1 atoms cm'2. The SN apparently lies beyond most of the neutral 
interstellar gas in the LMC (de Boer, Richtler and Savage 1987). 



From the information above we conclude it is likely that relatively strong to very strong A1 in, 
Si rV an C rV absorption will be produced by the pre-existing ISM in the direction of the SN. 
However, it is difficult to make a reliable estimate of the expected strength of this absorption 
because the ISM in the direction of the SN is so patchy. Some very rough guidance might be 
provided by lUE interstellar absorption line measurements toward the three objects R136a 
(de Boer, Koomneef and Savage 1980; de Boer, Fit 2 patrick and Savage 1985) , R144 (de Boer 
and Savage 1980), and HD36402 (deBoer and Nash 1982) whose A1 in. Si IV and C IV profiles 
are shown Figure 9. R136a lies at the center of the 30 Doradus nebula and has E(B-V)lj^C 
about 0.30. R 144 lies about 4.3 arc min from center of the nebula and has E(B-V)lmc about 

0.09 (Fitzpatrick and Savage 1984) . HD 36402 lies in a roughly spherical H H region of 
approximately 100 pc diameter located about 2° from R136a which places it well beyond the 
boundaries of the 30 Doradus nebula and has a small E(B-V)lj 4 C (see deBoer and Nash 1982) . In 
the progression from R136a to R144 to HD 36402, the Si IV and C IV absorption lines at LMC 
velocities progress from exceptionally strong to moderately strong. One might expect the data for 
the SN sight line to represent an intermediate case. However, we are unable to predict whether the 
absorption line strengths will lie closer to the exceedingly strong R136a case or the moderately 
strong R144 case. 

The expected velocity of the A1 III, Si IV and C IV absorption produced by H II region gas will 
depend on the gas kinematics near the site of the SN. It is well documented that the gas motions in 
the 30 Doradus region are violent and complex and there can be large variations from one position 
to the next (Elliott et al.l977; Meabum 1984). Y. H. Chu (private communication) has obtained a 
long slit high resolution echelle spectrum of a position 10" south of SN 1987A with an east-west 
slit orientation extending approximately 65" west and 130" east . The spectrum which includes the 
nebular emission lines of Ha X 6563 and [N II] X^6584 and 6548 was obtained 7 January 1988 
(UT) with the echelle spectrograph on the CTIO 4 m telescope, the air Schmidt camera and the 
GEC 385x576 element CCD. The strong telluric Ha emission feature gives us a zero point in the 
velocity scale and indicates that the 1.64" slit width provided a spectral resolution of 21 km s'^ 
FWHM. The SN measurements were part of a program to study SN remnants in the Magellanic 
Clouds. Details about the instrument configuration and the reduction techniques are given in Chu 
and Kennicutt (1988). The resulting spectrum shows strong Ha and weak [N II] emission from 
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d) Circumstellar Environment of the Supernova 

The immediate suiroundings of the supernova contain gaseous material expelled by the precursor 
star at vanous stages of evolution. When the star was a red supergiant, a slow wind probably created a 
large (-10 pc diameter) spherical volume of material centered on the star. Somewhere within the 
innermost parsec of this sphere, the much faster and lower density wind from the more recent blue 
supergiant phase swept up the wind from the cool star and compressed it into a thin, dense shell 
(Chevalier 1988). When the supernova exploded, the gas within this shell was ionized abruptly by the 
enormous impulse of ultraviolet photons created when the outermost portion of the supernova's 
envelope was heated to T ~ 1Q5.35 k by the emergence of the shock from the interior (Dopita et al. 
1987; Woosley 1988). Several months later, UV and visible emission lines from atoms in the densest 
portions of this compact shell became evident (Wamsteker et. al. 1987; Kirshner et. al. 1987) at a 
heliocentric velocity of approximately 280 km s' 1 and with a FWHM of about 30 km s' 1 (Fransson et 
al. 1988). The analysis of these data showed that the ejected material is especially rich in nitrogen - an 
indication that significant loss of the outer hydrogen envelope may have occurred in the star's earlier 
stages of evolution. Moreover, ratios of different emission line strengths indicated a characteristic 
density of electrons n^ = 2xl04 cm'3 (Fransson et al. 1988) and temperature T = 5.5x1 0^ K (Wampler 
and Richichi 1988). The angular extent of the region emitting the UV lines was not resolved by the lUE 
spectrograph, which indicates the radius of the dense shell must be less than 0.3 pc. The fact that some 
of the lines were still growing in strength about a year after the explosion (due to hght travel time 
effects) demonstrates that the shell must not be significantly smaller than 0.3 pc. Fransson et. al. gave a 
veiy approximate estimate for the mass of the emitting material; their value of 0.05 Mq is roughly 

consistent with the expected amount of gas from both the blue and red star winds within several tenths of 
a pc from the supernova. 

In principle, UV absoiption lines could provide addidonal information on the abundances and 
physical character of gases very near the supernova. Unfortunately, nearly all of these circumstellar 
lines are completely swamped by the very broad features created by the much larger amount of 
foreground gas associated with the LMC. Our only hope in bypassing the effects of the more pervasive 
gas at approximately the same doppler shift is to rely on absorption features from excited fine structure 
levels of vanous atoms. Under the right circumstances these absorption features may strongly favor the 
circumstellar gas, since the excited levels are populated either by collisions in regions of high density 
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and tcmperamre, or by optical pumping through higher electronic states in the presence of intense 

ambient radiation ftelds (Spitzer and Jenkins 1975). 

The most conspicuous UV lines from fine-stmctuie levels are those from the upper level of C n. In 
an H n region at a temperature 103 < T < Kk* K, an electron density ne - 30 cm‘3 is sufficient to create 
an excited C n population equal to that of the ground level. The relative ease in creating C n*. together 
»ith the high cosmic abundance of carbon, explains why we see C n* over a broad range of velociUes 
(more than 100 km s' 1 ) centered on 270 km s'l . Evidendy we are stUl registering the presence of 
foreground ionized (and possibly even neutral) gases associated with the LMC. 

Lines from Si II* are more difficult to produce. Within any H n region of low density, collisions 
with electrons will create a ratio of Si D* to Si U which is 17 times lower than that of C n* to C n -- a 
fact supported by the profiles shown in Figure 1. For optical pumping the disadvantage is much less. 

Si II has a number of very strong UV transitions from the ground state. It follows that Si D* should 
allow us to discriminate better the circumstellar absorpdons from features produced by the more 
extensively distributed material in the LMC. Hence, we shall rely on the Si if absorptions (or lack 
thereof) to probe the environment near the supernova. For our interpretation of the Si D* Unes, we first 
need to predict how many of the nearby Si atoms escape being ionized to higher levels by the U V 
photons from the supernova. A B3 la star produces l(f 5.57 Lyman limit photons s'l (Panagia 1973). 
This rate is about a factor 25 below that needed to maintain full ionization of 0.05 Mo's of hydrogen at a 
density ne = 2x1(7* cm'3. Therefore, provided the dense shell surrounding the supernova is reasonably 
intact (i.e., it has not broken into thin filaments with large gaps in between), most of the hydrogen 
within and outside the shell must have been neutral prior to the outburst 

Nearly all of the Si atoms coexisting with the neutral hydrogen would be singly ionized. Initially, 
the Si+ would be shielded from most of the supernova’s ionizing radiation (E > 16.34 eV) by 
intervening neutral hydrogen. At a distance r from the supernova, we expect that only after a time delay 

t(r) satisfying 


t(r) 


471 J J Ng (t) r J (t) dE dt = N„(r) 

0 13.60 eV 

will the ionization front eat its way through an interior volume contairang Nnfr) hydrogen atoms, e 
within the sheU or in the suirounding, undisturbed medium, and thereafter expose the Si II atoms to the 
ionizing photons. For simpUcity. we assume that the expanding supernova envelope radiates as a perfect 
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black body, yielding a photon energy spectrum 

QRQxlO^^E^ -3-1 -1. 

jsj/'t'i = (photons cm s erg ) 

exp[11605E/T(t)]-l 

To define the evolution of this flux with time, we adopted the effective temperature Tg(t) from 
Woosley's (1988) model lOH and calculated the envelope’s radius r* (t) = [L(t) / 47taTg^(t) 

Another simplification of our analysis is that the thickness of the front is assumed to be infinitesimal and 
that the hydrogen ahead of it is perfectly opaque to all photons with E > 13.60 eV. For the short time 
scale which applies to this problem, there are virtually no recombinations of the ions or hydrodynamic 
adjustments behind the front. 

Once the ionization front has past beyond any given point, the local Si"*" is exposed to all of the 
subsequent ionizing photons and has a survival probability, 

oo oo 

P(r) = exp 1^ - (r*/r)^ J J Ng(t) a(E) dE dt J (5) 

t(r) 16.34 eV 

where the photoionization cross sections a(E) for Si"*" are those calculated by Reilman and 
Manson(1979). 

Virtually all of the Si atoms within the dense shell are ionized by the LTV photons from the 
supernova. Beyond the shell, the undisturbed wind created long ago by the red supergiant has a 
hydrogen density 

n(H)=3.20xl0® (cm'^) , (6) 

2 

r V 

w 

where the mass loss rate M is expressed in M0 yr^ , r in pc, and the velocity of the wind v^^, in km 
. For any reasonable values of M and v^y, n(H) is too low to collisionally excite the Si II fine 
structure levels for r greater than a few tenths of a pc. However radiation from the supernova will excite 
the very strong ultraviolet transitions of Si II and optically pump the fine-structure levels of the ground 
electronic state. In the context that the pumping lines are optically thin, Flannery, Rybicki and Sarazin 
(1980) calculated that at a distance r from a spherical source of radius r* emitting Fy ergs cm'2 s'^ Hz'^ 
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an equilibrium ratio 

N(Si II*) ^ 0.678 % . (7) 

N(Si J 1 + 1.017 ^ 

will be established, where 

^ = 8.79x10'' . (8) 

(If the Si n is optically thick, the optically thin approximation must stiU apply in the inner part of the 
zone, and the thickness of this zone is sufficient to produce observable Si D lines if the velocity 
dispersion is not markedly less than the instrumental resolution of lUE.) The time constant for 
establishing the equilibrium is governed by the inverse of the spontaneous decay rate A 21 of the upper 
fine-structure level, which for Si U amounts to 1.3 hours. Since this time is short compared to the 
interval over which the flux could appreciably change at the time of observation, the equilibrium 
calculation should be valid. 

Near the middle of the time interval over which the shortest wavelength portions of the spectrum 
were observed (the first four images Usted in Table 1), Kirshner et al. (1987) found that the apparent 
flux of the supernova was about 1.1x10-23 erg cm'2 s'l Hz’l at wavelengths covering some of the 
strongest Si D lines. To convert this flux into an estimate at a representative distance of 1 pc, we need to 
compensate for the attenuation of the UV flux by foreground dust, in addition to multiplying the apparent 
flux by the square of the distance to the supernova, (5x10^ pc)2. Fitzpatrick (1988) esnmates that the 
total reddening E(B - V) toward the supernova is 0.15, of which 0.07 is created by material in our 
galaxy. If we assume E(1250A- V)/E(B - V) = 6.6 in our galaxy (Savage and Mathis 1979) and 8.0 m 
the LMCl (Fitzpatrick 1986) and consider that Av / E(B-V) = 3.1 in both regions, we calculate that 
"dereddening" the flux should amount to an increase by a factor 4.2, making the real flux at 1 pc equal to 
1.2x10-13 erg cm-2 s-1 Hz-1 and ^ = 1.0. 

The abundance A(Si) of atomic silicon to hydrogen in the circumstellar gas is probably considerably 
lower than ihe solar abundance ratio, since the ratio of heavy eletnenB to hydrogen in the LMC is lower 
than that of our galaxy. Also, some of Ihe Si may be depleted onto grains. Some insight on A(Si) in 


1 The value could be as large as 10 if the supernova has dust which is more 
of stars associated with the 30 Dor nebula. 


characteristic of that in front 
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LMC gases is available from our plot of the apparent optical depth of the line at 1808 A. An integration 
of X over the range 190 < v < 320 km s'l ( see Fig. 2). leads to a formal result for the column density 
N(Si II) = 3.510.3x10^5 cm*2. Unlike the profiles of the highly ionized species, however, w'e have no 
assurance that unresolved, deeply saturated structures within the profiles are not misleading us into 
underestimating the amount present. Moreover, many of the Si atoms in H II regions which are excited 
by very hot stars may be doubly or triply ionized. Thus, our inferred abundance A(Si) taken from 
N(Si n)/N(H) = 2x10'6 is actually a lower limit. Here, N(H) was taken to be 2xl021 cm‘2 (see §IVc). 

One might question whether or not A(Si) in the circumstellar wind could be different from the 
generally distributed material in the LMC because grains are created efficiently but the material has not 
yet been exposed to the usual destructive processes. On the one hand, a study of the mass-loss 
envelopes of a diverse collection of M giants in our galaxy by Knapp (1985) indicates that there are no 
anomalies in the gas-to-dust ratios. On the other hand, observations of very weak thermal SiO emission 
in the outer portions of such envelopes by Morris et al. (1979) indicated that perhaps less 
than 1% of the Si may have escaped being consolidated into grains. It is difficult to say whether or not 
this high efficiency in depleting the gas-phase Si is maintained for the modified abundances in the 
circumstellar gas of the supernova progenitor. 

Incorporating the principles discussed above, we arrive at our net expectation for the column density 
of Si n* through the immediate neighborhood of the supernova. 


N(SiII*) = 


J n(H)A(Si)P(r) 

r* 


N(Si n*) 


dr . 


(9) 


The greatest concentration of Si n* is at about 1 or 3 pc from the supernova; very few of the Si atoms 
remain in the singly ionized form within the dense shell. Table 5 summarizes the results of Eq. 9 for 
several values of MAw The numbers clearly show that the dependence of N(Si II*) on M / v^ is 
much steeper than linear. The principal cause of this effect is that the additional material slow's down the 
ionization front in the critical region where optical pumping is important. The resulting delay shields the 
Si atoms from the largest dose of hard, ionizing radiation which occurs at the very beginning. A factor 
of 2 change in ^ for appropriate values of M / v^ seems to change the resultant N(Si II*) by a relatively 
small amount. 
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The tracings in Fig. Id reveal a moderately significant absorption feature due to Si U* at 1264.74 A 
at a velocity near that of the material in the LMC. Corresponding features are not evident for the two 
other Si H* lines (KX 1309.27 and 1533.45 ) because these transitions are not as strong. Assuming that 
the 1264.74 A line is completely unsaturated, our measured equivalent width of 67 ±20 mA yields a 
column density N(Si II*) = 5.5xl0l2 cm'^ . The velocity centroid of the profile is at 255 ±6 km s'l. 
While we expect outflowing circumstellar material to have an absorption velocity slightly more negative 
than most of the gas which can be seen in emission, our displacement of -25 km s'l seems rather large. 

It is likely that nearly all of the Si H* we have registered belongs to a foreground H n region m the 
LMC, and indeed the velocity we measured is approximately consistent with this interpretation. Also, 
we calculate that at a temperature of 104 K, an H D region containing the amount of Si D* that we 
derived should have an emission measure of 600 cm’6 pc, assuming , as before, that SiAI = 2x10-6 ^ 
the LMC. The telluric Ha recorded by Y.H. Chu and shown in Fig. 10 can be used to extimate the 
intensity of the much stronger nebular Ha emission near 270 km s' A At the time of observation the 
expected intensity of the telluric line based on the Ha models of Ander«>n et al. (1987) is 
approximately 6.5 Rayleighs. This number assumes solar maximum conditions, which according to J. 
Harlander (private communication) provide better agreement between the model predictions and recent 
direct measures of the telluric Ha emission . Therefore, the nebular 270 km s*l Ha emission, which is 
25 times stronger than the telluric line, has an intensity corresponding to about 160 Rayleighs. 
Assuming that the Ha emission is optically thin and gives 0.243 Rayleighs for 1 cm 6 pc at T - 10 K 
(Osterbrock 1974) , this implies a nebular emission measure of 660 cm-6 pc. 

From the above argument, we propose that Si n*absorption from gas very near the supernova is of 
order or less than about half the total column density we detected at 255 km s" A For ^ = 1 , the 
numbers in Table 5 indicate that a conservative upper limit for any mass-loss material which is outside 
the dense shell gives a density less than that of a wind having M / v^ -lxlO‘6 m© yr 1 (km s A A 
Our upper limit rules out the upper end of the range of observed mass loss rates for cool, giant stars 
(e.g. Bowers 1985; Knapp 1985) such as those found for the OH/IR maser enussion stars (Herman 

1985). 

If the dense shell of gas which is responsible for the emission lines contains 0.05 M© at a distance 
of 0.3 pc, its column density of hydrogen N(H) would be 6xl0l8 cm'2 . On top of it one could add 
3 x 1019 cm-2 for the undisturbed wind material . Note that for our limit for M / v^ . the mass of 
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material piled up at r = 0.3 pc works out to 0.3Mq. From Eq. 6, the volume density just outside the 
shell should be < 35 cm*3 , with an inverse-square law decrease at greater distances. This amount of 
gas is substantially lower than the uniform value of 1000 cm-3 assumed by Raga (1987), who proposed 
that most of the C IV and Si FV appearing at the LMC velocity in the supernova spectrum was created by 
flash photoionization of the surrounding material. Quite apart from any specific model for the creation of 
Si rV very near the supernova, we find that if we multiply our infeired total column density for 
circumstellar gas by our adopted A(Si), we arrive at N(Si) = 7.2xl0l3 cm'2 - a value which is lower 
than our observed lower limit for N(Si IV) for gas at velocities near that of the supernova (see Table 3). 
Likewise, if the abundance of carbon in the circumstellar material is less than 1/10 the solar value, more 
than 40% of the C would need to be triply ionized if N(C IV) equaled our lower limit of lOl'^-^ ^ cm"2 
deduced from the 1550.77 A line. This statement is probably much too conservative, since the real 
column density of C IV is very likely to be substantially greater than our formal lower limit. 

In addition to absorptions from Si II*, we could in principle also look for evidence of excited 
neutral oxygen atoms. O I has two excited fine-structure levels, however for the 1302 A multiplet only 
the upper member of the two levels is free from interference by another line (the 1304 A line of Si II). 
The UV transitions which can pump O I are considerably weaker than those of Si II, however this is 
compensated by the much lower spontaneous decay rate A 21 . For a given amount of gaseous material, 
the visibility of O I is roughly equivalent to that of Si II, and the higher cosmic abundance of oxygen is 
balanced by the fact that the transition we can view is much weaker. The main drawback which makes 
O I** less sensitive than Si E* in our search for circumstellar material is that the statistical weight of the 
upper level of O I is only 1/9 of the total for all three levels, as opposed to 2/3 for the excited level of 
Si E. Thus, we do not gain by including O I** . 

Summarizing, our analysis of the essential features in the interstellar spectra at LMC velocities 
suggest that the absorption is dominated by gas which was present in the vicinity of Sk -69 202 before 
the star became SN 1979A. The strengths and the velocity structure in the lines of C IV, Si FV and Si E* 
are compatible with absorption by normal H E region gas. 


We thank the lUE observatory staffs at the Goddard Space Flight Center and at the ESA FUE ground 
station for their help in acquiring and processing the lUE satellite data for SN 1987A. The very rapid 
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estimating the expected strength of the terrestrial Ha line, a result we used in determining the nebular 
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TABLE 1 


HIGH DISPERSION lUE ECHELLE SPECTRA OF SN 1987A 


Image No. Exp. lime 
(min ) 

Date (UT lime) Ground^ 
(hr;min) Station 

wavelengths (velocity shift)'’ 
Angstroms (kms"') 

comment^ 

SWP 30377 

14 

Feb 24 (21:37) 

G 

1230 to 1800 (-4 );1800 to 1900 (0) 

image saturation for X > 1900 A 

SWP 30379 

14 

Feb 25 (00:53) 

G 

1230 to 1860(0) 

image saturation for X > 1860 A 

SWP 30381 

30 

Feb 25 (04:24) 

V 

1230 to 1800(0) 

image saturation for X > 1 800 A 

SWP 30383 

45 

Feb 25 (08:09) 

V 

1230 to 1800 (0) 

image saturation for X > 1800 A 

SWP 30384 

11 

Feb 25 (09:30) 

V 

1808 (-5); 1854 (-4); 1862(0) 

weak signal for X < 1700 A 

SWP 30389 

25 

Feb 25 (21:17) 

G 

1808 (-5): 1854 (-3) 

weak signal for X < 1700 A 

SWP 30394 

40 

Feb 26 (04:14) 

V 

1808 (-1); 1854 (+3); 1862(+5) 

weak signal for X < 1700 A 

SWP 30399 

75 

Feb 26 (13:44) 

V 

1808 (-4): 1854 (+5); 1862(+7) 

weak signal for X < 1700 A 

LWP 10190 

1.83 

Feb 24 (22:14) 

G 

2000 to 3100 (-18) 


LWP 10192 

1.67 

Feb 25 (00:45) 

G 

2000 to 3100 (-18) 


LWP 10194 

5 

Feb 25 (05:30) 

V 

2000 to 2700 (-10) 

image saturation for X > 2700 A 

LWP 10196 

8 

Feb 25 (07:52) 

V 

2000 to 2700 (-8) 

image saturation for X > 2700 A 


a Images from both the Goddard (G) and Vilspa (V) ground stations were utilized in this study. 


b We list the wavelength regions for which the various spectra were useful in our analysis and in parenthesis list the ty 
shifts in km s*^ needed to bring the various spectra into proper velocity registration . 


c TT^ diffaent spectra vary censideraWy in qnality and wavelength coverage. The SN faded very rapidly over the first few days 
of observations and reasonably complete wavelength coverage was only obtained with the early spectra. 


TABLE 2 


UV AbsoT>tion Lines Illustrated in Figure 1 


ion 

fa 

wavelength^ 

Echelle Figure 
order 

Notes 

Aim 

0.268 

1862.795 

74 

la 


Aim 

0.539 

1854.720 

74 

la 


Si IV 

0.262 

1402.769 

98 

la 


Si IV 

0.528 

1393.755 

99 

la 


CIV 

0.097 

1550.774 

89 

la 

1 

CIV 

0.194 

1548.202 

89 

la 

1 

OI 

0.0486 

1302.169 

106 

lb 

2 

Si II 

0.0055C 

1808.003 

76 

lb 


Sill 

0.147C 

1304.369 

106 

lb 

2 

Si II 

0.23C 

1526.719 

90 

lb 


Sill 

0.959 

1260.418 

109 

lb 


CII 

0.118 

1334.532 

103 

lb 

3 

Mg I 

1.77 

2852.127 

81 

Ic 


Mg II 

0.295 

2802.704 

82 

Ic 


Mg II 

0.592 

2795.528 

83 

Ic 

4 

AlII 

1.88 

1670.786 

82 

Ic 


Fen 

0.0395 

2373.733 

97 

Ic 


Fell 

0.203 

2599.395 

89 

Ic 

5 

Sin* 

0.860 

1264.730 

109 

Id 


Si n* 

0.147 

1309.274 

105 

Id 


Si II* 

0.076 

1533.445 

90 

Id 


Cl 

0.081 

1560.310 

88 

Id 

6 

Cl 

0.136 

1656.928 

83 

Id 

6 

Kill 

0.068 

1741.560 

79 

Id 



3 f values are from Morton and Smith (1973) except as noted, 
b X’s are in vacuum for k < 2000 and in air for X >2000. 


c f values from Shull, Snow and York (1981). 


notes to follow table 2 

NOTES... (,) AbsoT,Uon by boa, con,po„en.s of C IV double, appear op 
each plop Nop, d.a. a,e rligh, difference ip appearance in d« san.e 
componen. fron, one pio, ,o d,e „e„ ir d,e reap,, of binning u» daia a. 
slightly different velocities in each plot 

(2) Si n ^,304 absoipdon is ai,„aKd near d,e O I W302 abKnpPon. 

OIAbaorpdonby CII' W335.71 overlaps absorpPonbyCnW334.532. 

However. U« LMC C n* absonnion is free of eonp^inaPon. 

(4) Absoipuon by bod, lines of d,e Mg II doable, appear on each plot. 

(5) Tbefeamres.od.eIef,andrigb.ofFeII W599areMn„ 1,2593.73, and 

A2605.697, respectively. 


(6) The bnes of C I are acnially complex overlaping blends of absorpPon 

from C, .c- and c- See Morion and Smid. (1,73, for a bsPngof ab 
the components. 


TABLES 

COLUMN DENSITIES TOWARD ‘SN 1987A 


ION 

X(A) 

Interval 1 
( v<120km s"^) 
logN 

1.1. best^ u.l. 

Interval 2 

( 120<v<190 km s"^) 
logN 

1.1. best® u.l. 

— K 

Interval 3 
(v> 190km s*^ ) 
logN 

1.1. best® u.l. 

Aim 

1854.72 

13.09 13.12 13.15 

11.49 11.64 11.73 

13.45 13.49 13.53 

Aim 

1862.79 

13.11 13.15 13.18 

12.01 12.08 12.15 

13.48 13.51 13.54 

Si IV 

1393.76 

13.54 13.58 13.62 

12.72 12.77 12.82 

14.02 

Si IV 

1402.77 

13.48 13.52 13.57 

12.72 12.78 12.84 

14.19 

CIV 

1548.20 

14.10 14.15 14.21 

13.36 13.41 13.46 

14.59 

crv 

1550.77 

14.10 14.15 14.20 

13.40 13.46 13.51 

14.81 


a Lower limits 0-1.) . best values, and upper limits (u.l.) to column densiUes for each high ionization absorption line are 
listed for the three indicated velocity ranges. These estimates are based on a direct integration of the opucal depth profiles 
shown in Figures 2. 3 and 4. For the highest velocity range ( v >190 km s*l) . we only list lower limits for the Si IV and 

C IV column densities. Near v « 270 to 280 km s-1 the Si IV and C IV optical depths are so large they can not be reliably 
measured (see §IIIb). 



TABLE 4 


COLUMN DENSmES FOR HIGHLY lONlZED MILKY WAY HALO GAS 


Ion 

Predicted ^ 
NIsin bl 

Halo Stars^ 
< N Isin bl> 

HD5980C 
N Isin bl 

HD36402d 
N Isin bl 

SN 1987AC 
N Isin bl 

Aim 
Si IV 

crv 

NV 

OVI 

(3.3-6.4)x1012 

(4.3-7.9)x1013 

(2.8-3.6)xl0^^ 

(5.8-6.0)x1014 

«2x1013 

=lxl0l4 

=3x10^3 

>3x1013 

4.5x10^2 

2.8x10^3 

1.8x1014 

2.2x1013 

4.6x10^2 
1.8x1013 
5.4x1013 
< 1.8x10^3 

7.3x1012 

1.9x1013 

7.7x1013 


a Predicted column densities of highly ionized gas through the halo (N Sin Ibl ) on one side of the galaxy 
based on the time-dependent ionization calculations of Edgar and Chevalier (1986). The values assume a 
fountain mass flow rate of 4M(yyr on each side of the galactic plane. The mass flow rate has been 
adjusted to provide approximate agreement for the average observed strength of N V absorption toward 
high Izl halo stars. The two values of column density listed for each ion involve different assumptions 
about the sizes of the cooling regions. These same calculations are able to explain the 

C IV and O ID] emission measurements of Martin and Bowyer(1987) provided the fountain mass flow 
rate is about 2x larger than assumed above. 

*> Obsavations are from Savage and Massa (1987) for Si IV, C IV and N V and from Jenkins (1978) for 
O VI. The value for O VI is listed as a lower limit because the measures only extend to Izl * 1 kpc. 

The values listed are for the sight line to HD 5980 in the SMC from Fitzpatrick and Savage (1983) and 
refer to gas with v < 100 km s' 1 . 

d nie values listed are for the sight line to HD 36402 in the LMC from Savage and de Boer (1981) and 
refer to gas with v < 1 20 km s‘ 1 . 


e The listed column densities are for gas in the direction of SN 


1987A with v< 120 km s’ 1 (see Table 3). 



TABLES 


expected colum n densities of Si n* (cm*2) 

M / Vyv [ Mq yr"l (km 

5x10-8(4 = 2) 5x10-7(4 = 2) 1x10-6(4=1) 1x10 -6(4^2) 2x10-6(4=1) 2x10-6(4 = 2) 

4.9x1010 1.7x1012 5.6x1012 7.8xl0l2 j.OxlOls' 3.7xlOl3 



FIGURE CAPTIONS 

Fig. 1- Averaged high-resolution intd^tcllaiUncprofites for SN 1987A based on the lUE 
spectra listed in Table 1. In most cases, four spectra have been averaged to produce the prof.les 

shosvn. lUE flux has been plotted against heliocentric velocity for the lines listed in Table 2. The 
aero level of flux is indicated by the dc marks on the vertical axis. In some cases more than one 
ion absorbs in the special region illusnated (see the comments to Table 2). R indicates a detector 

reseaux. 

Fig. 2.- Plots of apparent log N(v) vs heliocentric velocity for three absorption lines of St fl: (1). 
M526.72 with log (K = 2.545 (solid circles), (2) 21304.37 with log SX = 2.283 (solid squares) , 
and (3) 21808.00 with log 12 = 2.545 (open circles). N(v) was defined according to Eqns. 1 and 
2, and is clearly not vahd for the two stronger lines, as explained in the text. Logarithms of the 
apparent optical depths T for the three lines are indicated by the scales on the right-hand side. 

Small symbols above and below the main ones indicate the range of possible systematic errors 
caused by uncertaintities in both the adopted continua and zem-intensity baselines (they do no, 
show errors attributable to each point's uncertainty caused by noise). 

Fig. 3.- Logarithmic plot of apparent optical depth T (as defined by eq. 1) vs hehocentric velocity 
for the A1 111 1 854.72 A line (solid circles) and 1 862.80 A hne (open circles). The points for the 
1862 A line have been moved upward by 0.3 to correct for the fact that its transition probability is 
smaller by a factor of 2. As with Fig. 2, small symbols on either side of the main points indicate 
the range of synemmic uncertainties. Log t = 0 on this diagram corresponds to N(v) = lOl 1.58 
cm-2 (km s-l)-l (seeEq. 2). 

Fig. 4.- Logarithmic plot of apparent optical depth x (as defined by eq. 1) vs heliocentric velocity 
for the Si IV 1393.76 Aline (solid circles) and 1862.80 A line (open circles). The points for the 
1402 A Une have been moved upward by 0.3 to correct for the fact that its transition probability is 
smaller by a factor of 2. As with Fig. 2, small symbols on either side of the main points Mcate 
the range of sysremnric uncertainties, log x = 0 on this diagram corresponds to N(v) = id l H 

cm"^ (km (see Eq. 2). 
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R*. 6.- log N(S. IV) - log N (A1 ID) versus heliocennic velocity derived fro ,u 
dep.hplo.sofFigures4and3 Re „ fr°n> the apparent optical 

gures and 3. Resulri are no. ploned for those ranges of velocity where the 
apparent optical depths are unreliable Si IV and AI m h 

innirir u . 1 V and AI in have very different line profiles whic 

pl.es fta. then colunu, density rados change subsBuhally wit 
figure. 
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with velocity as indicated in this 


F.g. 7.- log N(C IV) - log N(A1 DD versus hehoceuuic velocifr. derived fr u 

depd, ploB of Figures 5 and 3 Re, „ ved from the apparent optical 
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figure. 


g. 8. log N(C IV) - log N( Si IV) versus heliocentric velocity derived fr 
depd, ploB of Figures 5 and 4 R , 

* ’“''•‘'““'““^"°'P'P««<f<>r<l.oserangesofveloci.ywhereti,e 

apparent optical depths are unreliable. C IV and Si fV h • -i 

tf.eircoluo.od • '■“'"“'“'"“"Ppn.ftfes which imphes that 

tf.e.rcolumnde„s.D-ratiosarere,ativelyindepende„tofveloci.y. 



C Profiles for R136a, R 144, and HD 36402 from Savage and Meade 
(1989). Intensity is plotted versus heliocentric velocity with the zero level of intensity indicated by 
the tic marks on the vertical axis. R136a and R 144 are situated in the 30 Doradus nebula while 
HD 36402 lies beyond the boundaries of the nebula in its own H II region. The lines of A1 III, 

Si rV and C IV toward SN 1987A (see Fig. la) are intermediate in strength between what is seen 
toward R136a and R 144. The complex pattern of motions of the nebular gas in the 30 Doradus 
region and beyond probably explains the differences in velocity of the high ionization absorption 
from one sight line to another. When looking directly toward the center of the nebula (e.g. toward 
R136a) the outflow velocities are rather extreme and produce nebular absorption features which 
overlap the absorption near zero velocity due to the Milky Way. R indicates a detector reseaux and 
B indicates a bright pixel, most likely produced by a cosmic ray event. Each profile plotted 
represents an average of five spectra (for R136a and HD36402) or two spectra (for R 144). 

Fig. 10.- Hydrogen a X6562.808 emission line spectrum for a direction 10" south of SN 1987A. 
Intensity is plotted versus heliocentric velocity in km s'l. Weak telluric Ha is apparent along with 
the strong Ha emission produced by H II region gas in the general vicinity of the supernova. The 
continuum is probably light from the supernova scattered into the spectrograph This spectrum 
was obtained with the echelle spectrograph on CnO 4 m telescope by Y. H. Chu (private 
communication). 
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